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Tackling climate change is at the top of many agendas. In this context, emission trading schemes 
are considered as promising tools. The regulatory framework for an emission trading scheme in- 
troduces a market for emission allowances and creates a need for risk management by appropriate 
financial contracts. In this work, we address logical principles underlying their valuation. 
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1. Introduction 

The generic principle of an emission trading scheme is based on the so-called 'cap-and- 
trade' mechanism. In this framework, an authority allocates fully tradable credits among 
responsible institutions. At pre-settled compliance dates, each source must have enough 
allowances to cover all of its recorded emissions, or be subject to penalties. 

A mandatory cap-and-trade system involves its participants in a risky business with 
an obvious need for risk management. That is, certificate trading is usually accompanied 
by a secondary market for emission-related futures, including a rapidly growing variety 
of their derivatives. Their pricing is addressed in this work. 

Our contribution focuses on a methodology between equilibrium and risk-neutral ap- 
proaches. Due to the complexity of emissions markets, risk-neutral dynamics must be 
addressed in terms of explanatory variables, viewed as proxies of fundamental quanti- 
ties. Thus, we utilize equilibrium analysis to explain the role of fundamentals in risk- 
neutral allowance price formation. Thereby, the key issue is a feedback relation between 
allowance prices and abatement activity. Namely, we demonstrate that any increase in 
allowance price causes market participants to enforce emission saving in order to sell 
their allowances. Hence, an increasing allowance price encourages a supply of certificates 
and lowers the probability of non-compliance, which tends to bring down their prices. 
Apparently, the correct description of this feedback is the key to derivatives pricing. The 
present work focuses on this issue. On this account, our contribution goes beyond any 
risk-neutral approach to modeling of emission-related assets suggested in the existing 
literature to date. 
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2. Emissions markets 

The literature on this subject is enormous: it encompass hundreds of books and papers. 
For this reason, we focus only on those market models which are relevant in the present 
approach. 

Economic theory of allowance trading can be traced back to [8] and [14], whose au- 
thors proposed a market model for the public environmental goods described by tradable 
permits. 

Dynamic allowance trading is addressed in [7, 11, 13, 16, 17, 21, 22] and in the literature 
cited therein. 

Empirical evidence from existing markets is discussed in [9]. This paper suggests eco- 
nomic implications and hints at several ways to model spot and futures allowance prices, 
whose detailed interrelations are investigated in [23] and [24]. 



Econometric modeling is addressed in [1], where characteristic properties for financial 
time series are observed for prices of emission allowances from the mandatory European 
Scheme EU ETS. Furthermore, a Markov switch and AR-GARCH models are suggested. 
The work [15] also considers tail behavior and the hetcroscedastic dynamics in the returns 
of emission allowance prices. 

Dynamic price equilibrium and optimal market design are investigated in [2]. Based on 
this approach, [3] discusses the price formation for goods whose production is affected 
by emission regulations. In this setting, an equilibrium analysis confirms the existence 
of the so-called 'windfall profits' (sec [19]) and provides quantitative tools to analyze 
alternative market designs. 

Pricing of options was addressed only recently. The paper [6] discusses an endogenous 
emission permit price dynamics within equilibrium setting and elaborates on valuation 
of European option on emission allowances. The paper [18] and the dissertation [25] 
deal with the the risk-neutral allowance price formation within EU ETS. Here, utilizing 
equilibrium properties, the price evolution is treated in terms of marginal abatement 
costs and optimal stochastic control. Also, the work [5] is devoted to option pricing 
within EU ETS. The authors suppose that the drift of allowance spot prices is related to 
a hidden variable which describes the overall market position in allowance contracts and 
they make use of filtering techniques to derive option price formulas which reflect specific 
allowance banking regulations, valid in the EU ETS. Finally, the recent work [4] presents 
an approach where emission certificate futures are modeled in terms of deterministic time 
change applied to a certain class of interval- valued diffusion processes. 

The present work brings aspects of risk-aversion into the line of research followed in 
[3, 18] and [2], which we briefly sketch now. Within a stochastic model of an emissions 
market, a so-called central planer problem is introduced and discussed in [18]. Under 
additional assumptions, the authors formulate this problem in terms of continuous-time 
stochastic optimization. Furthermore, they provide economic arguments justifying why 
optimal control solutions correspond to an equilibrium of the emissions market. Inter- 
preting the allowance certificate price as the marginal abatement costs, particular explicit 
solutions are discussed and yield a dynamic stochastic model for allowance price evolu- 
tion. The work [2] starts from the opposite direction. In a discrete-time framework, the 
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Radner equilibrium of an emissions market is introduced and constructed via a solution 
of the central planer problem. The work [3] yields an extension: in a slightly different 
setting, it is proved that any market equilibrium is reached by this methodology. Thus, 
results from [2, 18] and [3] show that a quantitative analysis of emissions markets is 
tractable in terms of stochastic control theory. However, this connection is valid only if 
risk aversion is neglected, in other words, under the assumption that each agent possesses 
a linear utility function. Losing sight of risk aversion comes at the costs of unrealistic re- 
sults. Among other singularities, it turns out that the equilibrium allowance price follows 
a martingale (with respect to objective measure!) with the consequence that allowance 
trading can be arbitrary, only the final position must be adjusted accordingly. 

This work resolves all of these problems. Starting from the no-arbitrage property which 
is satisfied in an equilibrium of a market with risk-averse players, we show that the risk- 
ncutral allowance price dynamics exhibits the above feedback property, which we formal- 
ize as a fixed point equation, discussing its solution. Wc show that for such a risk-averse 
setting, our fixed point equation plays the same role as the central planer optimal control 
problem for the non-risk-averse situation. Namely, it provides a methodology to describe 
the market equilibrium in terms of aggregated quantities. However, this description is 
valid only from the viewpoint of the so-called risk neutral dynamics, not being suitable 
for discussing all interesting problems. Still, derivatives valuation is naturally addressed 
in, and can be obtained in, this setting. 

3. Mathematical model 

Let (f2, J 7 , P, (J-t)t=o) be a filtered probability space. Assume that is deterministic 
and agree that all processes considered in this work are adapted to (J r t )^_ . Write E*(-) 
and P f to denote, respectively, conditional expectation and conditional distribution with 
respect to Tt ■ Consider a market with a finite number I of the agents confronted with 
emission reduction. 

Emission dynamics. For each agent i S 7, introduce the stochastic process (El) t ~Q 
with the interpretation that E\ describes the total pollution of the agent i which is 
emitted within the time interval [t, t + 1] in the case of the so-called 'business-as-usual' 
scenario (where no abatement measure is applied) . Although each agent is considered as 
a potential producer, purely financial institutions are also covered with this approach by 
setting emissions to zero, that is, El = for t = 0, . . . , T — 1. 

Abatement. Consider the opportunity to reduce emissions. Each agent i can decide at 
any time t = 0, . . . , T — 1 to reduce its emissions within [t, t+1] by Q. pollutant units. We 
suppose that each abatement level is possible, ranging from no reduction to full reduction. 
Hence, we assume that < £J < El holds for all t = 0, . . . , T — 1. 

Abatement costs. We assume that the cost of abatement is a random function of the 
reduced volume. The randomness is due to uncertainty in prices (of fuel) and is observable 
at the corresponding time. Thus, if the agent i decides at time t = 0, . . . , T — 1 on reduction 
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of their own emissions by x € [0,oo[ units, then it causes costs C\{x) 1 where given 

C\ : [0, oo[xf2 1-)- R is B([Q, oo[) ® J" t -measurable 
and for each u> £ Q,x h- > Cl(x)(u)) is strictly (3-1) 
convex and continuous with C(0) = 0. 

Since emission savings cannot exceed the business-as-usual emission, the abatement ac- 
tivity (£.1)^=0 is feasible if 

0<Q<EI f = 0,...,T-l. (3.2) 

Following abatement policy (£l)Js , the agent i accumulates at the compliance date T 
the total terminal costs 

T-1 
t=0 

Abatement volume. For later use, let us introduce, for each we!], t = 0, . . . ,T — 1 and 
a G [0,oo [, the abatement volume cj(a)(w) as 

cj(a)(w) :=argmin{Q(x)(ij) -ax:x£ [0, E l t (uj)]}, (3.4) 

which is well defined since, under the assumptions (3.1), the minimum of the function 
x i — y CI(x)(uj) — ax on [0, EI(uj)] is attained at the unique point. The reader may imagine 
cj(a)(a;) as the total reduction volume which is available within [t,t+ 1] in the situation 
w at a price which is less than or equal to a (measured in currency unit per pollutant 
unit). A straightforward proof shows that (3.1) ensures that 



(3.5) 



[0, oo[h-> R, a h-> c\{a){uj) is non- decreasing and 
continuous for almost every uj Gfl and t = 0, . . . , T — 1. 

For later use, we introduce the cumulative abatement volume function 

c t (a):=Y, e t{a), oe[0,oo[. (3.6) 

Obviously, Ct(a)(u>) stands for the total abatement in the market, which is available from 
all measures in the situation uj whose price is less than or equal to a £ [0, oo[. 

Allowance trading. Suppose that, at any time t — 0, . . . ,T, credits can be exchanged 
between agents by trading at the spot price A t . Denote by "&1 the change at time t 
in allowance number held by agent i. That is, given the allowance prices (A t )J =0 , the 
position changes (i9J)^ yield costs 

T 

5>jA t . (3.7) 

t=o 
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Penalty payment. The total pollution of the agent i can be expressed as a difference 

T-1 T-1 
i=0 t=0 

of the cumulative business-as-usual emission less the entire reduction. As mentioned 
above, a penalty ir s]0, oo[ is being paid at maturity T for each unit of pollutant, which 
is not covered by allowances. Considering the total change in the allowance position 
J2t=o^t e ff ecte d by trading, the loss of the agent i resulting from potential penalty 
payment is 

/T-1 \ + 

^n2(Ei-e t -^)-Y~^Tj , (3.8) 

where 

j l , i € I are agents' initial allowance allocations. (3-9) 



Remark 1. Our stylized scheme deals with stand-alone emission trading mechanisms. 
In the real world, cap-and-trade systems operate on multi-period scales, where unused 
allowances can be carried out (banked) into next period. Further period interconnec- 
tions may include a transfer of future allocation from the next into the present period 
(borrowing) and, in the case of non-compliance, a withdrawal of an appropriate number 
of credits from the next period allocation in addition to penalty payment. To complete 
the complexity, let us mention that different emissions markets could be interconnected 
by acceptance of foreign certificates in the national scheme. Emission trading in multi- 
period settings is addressed in, among others, [4] and [5]. Mathematically, it reduces to 
the specification of a more complex penalty mechanism than that presented above. For 
this reason, we have decided to focus on the stand-alone allowance market to analyze 
quantitative methods in the simplest situation before tackling multi-scale systems (such 
as the second period of EU ETS). 

Recording uncertainty. In what follows, we also need to take into account uncertainty in 
the emission recording. It is convenient to subtract these recording errors from the initial 
allocation. Hence, we interpret Y as the credits allocated to the agent i less emissions 
which become known with certainty only at time T. With this interpretation, 7^ stands 
for allowances effectively available for compliance and is modeled by an J^-mcasurable 
random variable. For later use, let us agree that the distribution of conditioned 
on J-t-i, possesses almost surely no point masses, which implies that 

P[^2f = X j = for each J" T _i-mcasurablc X. (3.10) 
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Admissible policies. Since maximally possible reduction cannot exceed emission, we 
have (3.2). Let us define the space of feasible trading i} 1 = (i9|)|L and abatement strate- 
gies = (Q) t =o of the agent i G I by 

:={(^,f):0<g<^,f = 0,...,T-l}. (3.11) 

Individual wealth. In view of (3.3), (3.7) and (3.8), the revenue of the agent i following 
admissible policy ($%£*) equals 

T-l 



L A > i (# i ,e)=-j2(4A t +c i (&) 



fT-l 



, t=o 




(3.12) 



i?isfc aversion. To face risk preferences, suppose that attitudes of the agents i 6 / 
are described by utility functions U l : R i— ► R, which are continuous, strictly increasing 
and concave. Consider the utility functional u z (X) = E(U l (X)), which is assumed to 
be defined for each random variable X where the expectation is finite or +oo. Given 
allowance price process A = (A t )f =0 , the agent i behaves rationally, maximizing (d 1 , £ J ) i— > 
u % (L A '' l ('d' 1 , £ 1 )) by an appropriate choice of their own policy (^",^ ! *). 

Market equilibrium. Following standard theory, a realistic market state is described by 
the so-called equilibrium - a situation where the allowance price, positions and abatement 
measures are such that each agent is satisfied by their own policy and, at the same time, 
natural restrictions arc fulfilled. In our framework, an appropriate notion of equilibrium 
is given as follows. 

Definition 1. The process A* = {A^)J =0 is called an equilibrium allowance price process 
if, for each i G I, there exists G IA 1 such that u l (L A * ' l (fl* 1 , £**)) is finite and 

(i) the cumulative changes in positions are in zero net supply, that is, 

^0* = O for allt = 0,...,T; (3.13) 

(ii) each agent i El is satisfied by their own policy, in the sense that 

u l (L A *' l {r\C 1 )) > u l {L A "' l {d\C)) 



for each ) where u x (L A * > % {f )) exists 



(3.14) 



The existence of emissions market equilibrium is addressed in [2] and [3], under the 
assumption of a linear utility function and in a slightly different setting. However, al- 
though equilibrium modeling in the spirit of these contributions is appropriate to inves- 
tigate important questions of optimal market design, it has little to offer to the problem 
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of derivatives valuation. With the present approach, we intend to establish a reduced- 
form model which describes the evolution of emission-related assets from a risk-neutral 
perspective. We obtain a realistic picture by incorporating three essential assumptions 
into a risk-neutral model. These assumptions are shown to be direct consequences of an 
equilibrium situation: 

(a) There is no arbitrage since, in equilibrium, any profitable strategy would imme- 
diately be followed by all agents. This would instantaneously change prices and 
exhaust any arbitrage opportunity. 

(b) The allowance trading instantaneously triggers all abatement measures whose costs 
are below allowance price. The explanation here is that if an agent possess a 
technology with lower reduction costs than the present allowance price, then it is 
optimal for that agent to immediately reduce pollution and take profit from selling 
allowances. 

(c) There are only two final outcomes for allowance price. Either the terminal al- 
lowance price drops to zero or it approaches the penalty level. The reason is that 
at maturity, the price must vanish if there is an excess in allowances, whereas in 
the case of their shortage, the price will rise, reaching penalty. We believe that 
in reality, an exact coincidence of allowance demand and supply occurs with zero 
probability and can be neglected. 

Let us formalize the above assertions (a), (b) and (c). 

Proposition 1. Suppose that (Af)J =0 is an equilibrium allowance price and (Q*)JSq 
for i 6 I are the corresponding equilibrium abatement policies. 

(a) There exists a measure Q which is equivalent to P such that (A^)f =0 follows a 
Q -martingale. 

(b) For each i £ I , we have 

£f=<4(A*), i = 0,...,T-l, (3.15) 

with abatement volume functions c\, t = 0, ...,T — 1, from (3.4). 

(c) The terminal value of the allowance price is given by 

a t = ,rl {E 16 xCC£-o 1 (^-eJ*)-7 4 )>o}' ( 3 - 16 ) 

Before we proceed with the proof, let us emphasize that this result can serve as a 
starting point for risk-neutral modeling. The above proposition states that at equilibrium, 
the allowance price process (Al)J =0 follows a martingale with respect to an equivalent 
measure Q ~ P whose terminal value is 

^r = 7rl{£ ieJ £?=o(^-q')-7 i )>0}> 

obviously depending on intermediate values (A% through abatement volume function 
Q* = c\{A* t ) for t = 0, . . . ,T— l,i G I. The surprising and far-reaching consequence is that, 
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from a risk-neutral perspective, only cumulative market quantities are relevant. To see 
this, define the overall allowance shortage 




(3.17) 



which would appear in the market without any emission penalty. Further, recall from 
(3.4) and (3.6) the cumulative abatement functions to express the risk-neutral certificate 
price dynamics in terms of the following feedback equation: 



t = 0,...,T-l. 



Although individual market attributes and actions of the different agents seem to be 
irrelevant in this picture, the reader should notice that this picture appears only from 
the risk-neutral viewpoint. In line with standard aggregation theorems, the equilibrium 
market state heavily depends on, and is determined by, market architecture, rules, risk 
attitudes and uncertainty. However, once equilibrium is reached and all arbitrage op- 
portunities arc exhausted, asset dynamics can be considered under risk-neutral measure. 
With respect to this measure, market evolution appears as if it were driven by cumulative 
quantities only. 

With this in mind, let us formulate the problem of the reduced-form modeling as 
follows: 



Given measure Q <~ P, random variable Et 
and abatement volume functions {ct)t=o ■> 
determine a Q-martingale {A^)J ={) with 



(3.18) 



7Tl 



{St-T.JS, 1 c t (At)>0}- 



Note that this formulation serves as a guideline for martingale modeling since price- 
dependent abatement volume c f (a) can be estimated from market data, whereas potential 
allowance shortage &r can be modeled in terms of total allowance allocation and demand 
fluctuations on goods whose production causes the pollution. Finally, we shall emphasize 
a natural passage to continuous time. 



' Given, on a probability space (f2, T, P, (•7 r t)tG[o,T])i 
an equivalent measure Q ~ P, random variable Et 
and a family of abatement functions {ct)te[a,T\i 
determine a Q-martingale (A% )te[o,T] with 



(3.19) 



■ irl 



{£ T -J^c t (At) dis- 



proof of Proposition 1. (a) According to the first fundamental theorem of asset pricing 
(see [10]), it suffices to verify that if (Al)f =0 is an equilibrium allowance price process, 
then there is no arbitrage for allowance trading. Let us follow an indirect proof, supposing 
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that iytj^^Q is an allowance trading arbitrage, meaning that 



(T-\ \ /T-l 



P [Y^Vt{At+i-A t )>0 =1, P [J2^(A t+1 -A t )>0 >0. (3.20) 



\t=o / \t=o / 

Now, we verify that in the presence of arbitrage, no equilibrium can exist since each 
agent i can change their own policy , to an improved strategy satisfying 

u i (X A '' i (r,f i ))<« i (L j4 *' i (iJ i ! f i )). (3.21) 

The improvement is achieved by incorporating arbitrage {vt)t=Q into their own allowance 
trading as follows: 

$\ := + (v\ - v\_ x ) for alH = 0, . . . ,T, 

with appropriate definitions v-\ = vt '■= 0. Indeed, the revenue improvement from al- 
lowance trading is 

T T T-l 

t=0 t=0 i=0 

which we combine with (3.20) to see that there is no optimality since 

p(i' 4 ' ! (r I e 1 )<i A ' ! (i? , ,f)) = i, p(L A ' % {$*\c ) <l a ^{¥,c )) >o 

together imply that 

u i(L A <\r\e))<u i (L A < i {&,e)). 
(b) To prove (3.15), consider the bijection 

W^U\ (8\e)^ (#?\0,f), (3.22) 

where the transformed trading strategy d 1 = <p(6 l , £*) is given by 

0j = ^-& i = l,...,T-l, 4 = ^. 

Obviously, (t? 1 *,^ 1 *) is a maximizer to the original problem 

if and only if (tf 1 *^") = (0(0",£")jC*)j where (#",£") is a maximizer to the trans- 
formed problem 

U l ^R, (6\C)^u l (L A '*' l ((j)(9\C),C))- (3-23) 
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The last line in the calculation 

T-1 T-1 /T-1 \ + 

= - E( * - £q(^) - W E(^-s - («*-s))-7 < -^t 

t=0 t=0 \t=0 / 

T-1 /T-1 \ + T-1 

= -J2eiA* t -Trlj2(Ei-ei)-l i -ei r ) -A t *£) (3.24) 

t=o \ t=o ) t=o 

shows that if (#",£") is a maximizer to (3.23), then £* must satisfy Q* := c\{Al) for 
t = 0, ...,T— 1, which proves (3.15). 

(c) This assertion is proved by an argument identical to that given in [2]. □ 



4. Reduced-form modeling 

In what follows, we propose a solution to the problem of risk-neutral allowance price 
modeling (3.18). Below, we prove that under the assumptions given above ((3.10), in 
particular, is essential), the problem (3.18) possess a solution. Moreover, we show how 
to obtain the martingale (Af )f =0 . 

It turns out that the martingale closed by £t plays a crucial role, so we introduce 

£ t =E e *(£ T \J 7 t), t = Q,...,T. 

For later use, let us also define its increments as 

Et=£t— £t-ii t = l,...,T. 

Following the intuition that the equilibrium allowance price should be uniquely deter- 
mined by the present time and the general market situation, we express a candidate for 
allowance price as 

AZ(u)=at(G t (u)){w), weO, t = 0,...,T, (4.1) 
with hypothetic functionals 

a t :lxfi^[0,4 B(E)(8)J r t -measurable for t = 0,...T, (4.2) 

applied to 

t-i 

G t = £ t -J2c s (A* s ), t = 0,...,T. (4.3) 

According to (3.18), this approach yields an obvious definition for or: 

a T (g)(u) = irl[o, 00 [(g), ujen, g eR. (4.4) 
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Note that, given functional (4.2), the price process (A^)f =0 is indeed well defined by 



recursive application of (4.3) and (4.1): 

A* t (Lj):=a t (G t (u;))(Lj), (4.5) 

G t +i(w) := G t (u) - ct(AZ(u)){w)+e t +i(w), (4.6) 

started at G := £q. (4.7) 



Generated by this recursion, the process (Aj)^L follows a martingale if, for all t = 
0, . . .T — 1 and almost all lo G f2, the following holds: 

a t (<7)(w)=E t Q (a t+ i( 5 -c t (a t ( ff )(w))(o;) + £ t+ i))(a;) for all jet, weSl. 

Indeed, we have 

E?( J 4* +1 )H=E?(« t+1 (C7 t+1 ))( w ) 

= / a t+1 (C7 t ( W ')-c t (An^))(^')+£ t +i^'))(^)Q t (d^)H 
Jo 

= / a t+1 (G 4 (o;)-c t (A*(a;))(a;)+£ t+1 (wO)(a; , )Qt(dcj / )M 
isi 

(w) - c t (at(G t (w))(w))(a;) + £ t+ i))(w) 
= a t (Gt(w))(w)=^(a;). 

In other words, it is sufficient to ensure that 

for each g G R, at(g)(a;) solves 

a = E?(a t+ i{g-c{a)+e t+1 )){u>) (4.8) 
for almost all u G O. 

In the remainder of this section, we will show that the functionals (4.2) are recursively 
obtained as the unique solution to (4.8), starting with o.t from (4.4). First, let us prepare 
an auxiliary result dealing with the solution to (4.8) where no conditional information 



needs to be considered. 
Lemma 1. Given 

c:R— ^R, non- decreasing, continuous, (4-9) 

ax :R x Q, ->■ [0,tt], B(R) <£> J 7 -measurable, (4.10) 

g i — y a\(g)(oj), non- decreasing for almost all w G f2, (4-11) 
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suppose that the random variable £ satisfies 

R-»[0,7r], x^E Q (a 1 (x + e))= / ai(x + £(o/))(o/)Q(du/) 

Jn 

is continuous. 

For each j£R, introduce the function f 9 : [0, it] — » R given by 
F(a) :=a-E Q (ai{g-c{a)+e)) 

= a- I a 1 {g-c(a)+e{uj r ))(uj')q(duj'), ae[0,ir]. 
Jn 



(4.12) 



(4.13) 



The following assertions then hold: 

(i) for each j£R, i/iere exists a unique ao(g) € [0,7r] / 9 (oto(<7)) = 0; 

(ii) root ao(g) o/ / 9 is obtained as a limit ao(g) = Iim,,-^ m f/ie standard 
bisection method 

" 2Kl+ - J ' < +1 =a«, < +1 :=<C 0, l4 ' i4j 



started at Oq := 0. afj := 7r: 
(iii) £/ie mapping K — > [0,7r], j4 ao(<7) «s non- decreasing and continuous. 



Proof, (i) For each <? € R+, the function f 9 is continuous due to (4.12) and the continuity 
(4.9) of c. Thus, the existence of a root follows from the intermediate value theorem 
because of 

/ ff (0)<0, /»>0. (4.15) 

The uniqueness of the root is ensured by the strict monotonic increase of f 9 . To verify 
this, observe that (4.11) and (4.9) imply that the subtrahend 

ax(g - c(a) + e(w'))(w')Q(du/) 



Jn 

in (4.13) is non- increasing, whereas the minuend a n> a is strictly increasing in a. 

(ii) The bisection algorithm is properly initialized because of (4.15). Standard argu- 
ments ensure its convergence to the root. 

(iii) To show the monotonic increase of g <— > ao(g), suppose that g' < g. Then (4.11) 
ensures that for each a£ [0,7r], 

a 1 ( 5 '- C (a)+ £ ( w '))(w')Q(d^)< / ai(5-c(a)+e(w'))(^)Q(dwO: 
q Jn 

giving f 9 (a) > f 9 (a) for all a & [0,7r], which implies that ao(g') < ao(g). 
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Now, let us turn to the continuity. If ao(g) € [0, tt[, then there exists S > with ao(g) + S < 
tt. Due to the strict monotonic increase of f 9 , we obtain < f 9 (ao{g) + 6). If (g n ) n eN C 
R + is a sequence with lim„_s. 0O g„ = g, then according to (4.12), 

lim f^(a (g)+6)=f g (a (g)+6)>0. (4.16) 

71— >0O 

Hence, there exists TV 6 N such that f 9n (ao(g) + 5) > holds for all n> N. Thus, we 
obtain 

a (ff) e [0,7r[ =*> limsupa (c/„) < a (g) + 5, if a (g) + 6 <tt. (4.17) 

n— >oc 

Since 5 > is arbitrarily small and < cto(g) <tt, due to (i), this implication shows that 
ao(-) is continuous on each point g with ao(g) — 0. A similar argument yields 

ao(ff)e]0,7r] =*> liminfa (ff™) >a (5) ~<5, ifa (flr)-5>0. (4.18) 

Again, since 5 > is arbitrary, we obtain the continuity of «o(') on each point g with 
ao(g) = tt. If ao(g) G]0,7r[, then the continuity of ao(-) on g follows by the combination 
of (4.17) and (4.18). □ 

Let us now turn to the conditioned version of Lemma 1. Supposing the existence of 
the regular /(-conditioned distribution Q t , the proof reproduces the arguments of the 
previous lemma with appropriate notational changes due to conditioning on the event 
lo <E tt. However, a useful insight is that the approximating points a® , n = 0, 1,2, . . . , of 
the bisection algorithm turn out to be dependent on g G R and lu <G SI in a S(R) ® ,F t - 
measurable way, which shows that the functional under discussion, (g,u) H> ctt(g)(oj), is 
also S(R) ® J^-mcasurablc, being the limit of the sequence ((g,u>) i— > a^' w )^Lo of measur- 
able functions. 

Lemma 2. Suppose that for t S {0, . . . , T — 1} , 

c:Ix!l->I, S(R) ® F -measurable such that (4.19) 

a4C((a)(ki) is non- decreasing, continuous, (4.20) 

a t+ i :R x O -> [0,7r], B(R) ® T -measurable such that (4-21) 

g I— > at+i (5) (a;) is non- decreasing for all ajgSl. (4.22) 

Given a regular version Qt of the J- 1- conditioned distribution Q, assume that the random 
variable £t+\ satisfies 



R "►[(), tt], / a t+ i(:E + e t+ i(w'))(w')Qt(da/)(w) 

° (4-23) 
is continuous for each wgfl. 

TTie following assertions then hold: 
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(i) there exists a unique £>(R) ® J-'t -measurable [0,ir]-valued at satisfying 

at(g)(u)=E?(at + i(g-ct{at(g))+e t+1 )){u) 

(4.24) 

for all j£t, for almost all lo e f2; 

(ii) the mapping R — > [0,7r], j >4 af(j)(u) is non- decreasing and continuous for all 

uj eft. 



Proof, (i) As in the proof of the Lemma 1, wc obtain the unique root at(g)(uj) of the 
function 



f 9 > u (a):=a- / a t+1 (. 9 -ct(a)H+e t+1 ( W '))(^)Qt(dw')H, aG[0,7r]. 
Jo 

By the bisection method, 



g,cj _ 1/779," | „g,w\ n+1 n 



= o(<' + <' ), 



a n+l — a n' U i Q-ri+1 -~ a ni H f 9 ,LJ ) < , 



started at Og' := 0, ap'" := tt. Since 

(g,ui,a)^> f 9 > u (a) is 0(R) ® Ji ® £([0,7r])-measurable, 

each bisection point (g, w) n- a 9 ^ is 0(R) <£> J^-measurable, which shows that for n — > 
oo, the pointwise limit (<?,w) H> ctt(g,w) of the bisection sequence is also £>(R) ® J~t- 
mcasurable. By construction, the equality 

at(s)M= / "t+iCfl 1 -Cf(«t(5)M)(w) +£ t +i(w'))(w / )Q«(dw / )(w) 



holds for all g € R and a; € f2, whose right-hand side is nothing but the right-hand side 
of (4.24) for each jel. 

(ii) The proof is obtained from (iii) of the previous lemma by replacing ai(-), cto(-), 
c(-) and Q(du/) by at+i(-)(w),at(-)(ui), Ct(-)(u>) and Qt((ko')(u>), respectively, with ap- 
propriate notational adaptations according to the conditioning on u. □ 

Finally, we address a solution to (3.18) in the last point of the following proposition. 

Proposition 2. Consider Et under the model assumption (3.10) and the cumulative 
abatement volume functions from (3.6) under (3.1) and (3. 4). 

(i) Given measure Q ~ P, there exist functionals 

a 4 :R xfi ->■[(), tt], S(R) <g> T t -measurable for t = 0, . . .T, (4.25) 
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which fulfill, for all g 6K. 

Ot(s) = 7rl[o,oo[(ff)i ( 4 - 26 ) 
a t ( ff )=E t Q (a t+ i( 5 -c t (a t ( 5 ))+£ t+1 )) ) t = 0,...,T- 1. (4.27) 
(ii) There exists a Q-martingale (A%)J =0 which satisfies 

A T = 7rl {£ T -^r = - 1 c t (Ai)>0}- ( 4 - 28 ) 

Proof, (i) In this proof, we repeatedly make use of Lemma 2. Let us start with t = T—l 
and verify that the assumptions of this lemma are satisfied. Due to continuity (3.6) of the 
abatement function, we have (4.9). The properties (4.21) and (4.22) hold for t = T — 1, 
by definition (4.26). To show (4.24), wc utilize the specific form of ar- 



x^ / a T (x+e T (uj')){u')QT-i{duj'){cv)=QT-i(x + £T>0){u)- (4.29) 
Note that, due to (3.10), there are almost surely no point masses in the distribution of 



conditioned on Tt-i (with respect Q, since Q ~ P). That is, (4.29) is continuous for each 
ui S fi, as required in (4.24). Hence, (i) of Lemma 2 yields the functional a-r-i satisfying 
(4.27) (with t = T — 1), as required. To proceed by induction, we emphasize that (ii) of 
Lemma 2 ensures that g i— > aT-i(<?)(w) is non-decreasing and continuous for all to £ f2. 
That is, for the next step, i = T — 2, the assumption (4.22) on q:t-i is automatically 
satisfied. Moreover, (4.24) now follows, due to the continuity of g i-> ar-i(g)(w), from 
the pointwisc convergence 

lim ctT-i(x n + £T-i(^'))( a; ') — a T-i{x + £t_i(u/))(o/) for all a;' e 

n— >c« 

dominated by 7r, which holds for each (x n ) nS N C M with linin^ooXn = a;. That is, all 
assumptions of Lemma 2 are also fulfilled for t = T — 2. Proceeding recursively for t = 

T-2,...,0, we obtain (a t )^ with ( 4 - 25 ): ( 4 - 26 ) and ( 4 - 27 )- 
(ii) As suggested by (4.5)-(4.7), we define, for all we!!, 

A*(u) :=a t (G t H)H, 
Gt+i(u) := G t (u) - c t (A;(u))(w)+e t+1 (uj), 
started at Go := So- 

The process (y4|)^_ generated in this way obeys the terminal condition (4.28), in view of 
(4.26). To show the Q-martingale property of (A*)f_ , we calculate, for t = 0, ...,T— 1, 

E?(^ +1 )H=E?( at+1 (G t+1 ))H 
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= Ef(a t+1 (G t (uj) - c t (A* t (u,))(u>) + e t+1 ))(u,) 
= E?(a t+1 (G t (u) - c t (a t (G t (w))(iu))(uj)+e t+1 ))(oj) 
= a t (G t +i(u>)){uj) = A* t (uj) 
for almost all w £ f2, where the penultimate equality follows from (4.27). □ 

5. Applications 

Let us elaborate on the computational feasibility of our reduced-form modeling. For 
illustrative purposes, we focus on the simplest case of martingales with independent 
increments and deterministic abatement functions. We assume that: 

Et+i and Ft are independent under Q for all t = 0, . . . ,T — 1; (5-1) 
Ct : [0, oo[— } K is deterministic and time constant (ct = c)^ 1 . (5-2) 

Under these assumptions, the randomness enters the allowance price through the present 
up-to-day emissions only. More precisely, (5.1) ensures that 

uj i y at(g)(uj) = at{g) is constant on f2. (5-3) 

Let us verify this assertion. For t — T, (5.3) holds, by definition (4.26). For t = T— 1, . , . , 1, 
we proceed inductively as follows: by construction, a t (g)(uj) is the unique solution a to 

a= a t+1 (g-c t (a)+e t +i(Lu')){uj')Clt{du')(u) 
Jn 

= / a t+ i(.g-Q(a)+e t+ i(u/))Q(du/), (5.4) 

where, in the last equality, we have utilized the fact that Qt = Q, due to the independence 
(5.1), and the fact that at+i{g) does not depend on cj, by the the induction assumption. 
Obviously, the fixed point a t (g)(uj) := a from (5.4) also does not depend on uj. 

For numerical calculation, we rely on the one-dimensional least-squares Monte Carlo 
method, which is applicable in our case of martingales with independent increments. 
Although this setting is relatively restrictive, it covers a sufficiently rich class of mar- 
tingales. For instance, important cases of information shocks leading to allowance price 
jumps can be easily addressed under this approach when (£t)f = o is modeled as an ap- 
propriately sampled, centered Poisson process. In this case, fixed point equations can be 
treated analytically. We do not follow this path in favor of numerical methods, which 
deserve particular attention due to the complexity of emissions markets. In particular, 
extensions of Monte Carlo methods to the multidimensional setting (see [20]) seem to 
be appropriate. A preliminary analysis shows that assuming the existence of a global 
Markovian state process allows independence to be weakened to conditional indepen- 
dence, which leads to multidimensional Monte Carlo, in the sense of [20], since the state 
process gives additional dimensions. 
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We now focus on computational aspects. From (5.3), it follows that a t {G t ) is a <r(G t )- 
measurable random variable. Thus, in the equality (4.27), the condition J- t can be re- 
placed by the condition a(G t ): 

a t (G t )=EQ(a t+1 (G t -c t (a t (G t ))+e t+1 )\<j(G t )). (5.5) 

We shall treat this relation as a fixed point equation for the Borel-measurable function 
at and attempt to obtain a solution in the limit at = limn^oo a" of iterations 

a^ + \G t )=E Q {oL t+l {G t -ct{a^{G t ))+e t+l )\cx{G t )), neN, (5.6) 

started at a® = at+i- (Note that, given at+i and a", the equation (5.6) indeed defines 
a Borcl function a™ +1 by the factorization of the (j(Gi )-measurable random variable on 
the right-hand side of (5.6).) For numerical calculation of conditional expectations, we 
suggest using the least-squares Monte Carlo method. 

To explain the principle of the least-squares Monte Carlo approach (see [12] and [20]) 
in more detail, we abstract from the concrete situation (5.6) and consider 

^(G)=E°W(G,e)KG)), 

where e, G are K- valued and independent with respect to Q and <f> is a bounded 
Borel function on R 2 . Under these assumptions, the function ip is obtained as ip(g) = 
J K 4>(g,e)Q £ (de) for Q G -almost all g e R, where Q £ ,Q° are image measures of Q under 
e and G, respectively. An equivalent condition defining ip is the orthogonality 

determine ip £ L 2 (R, u) such that for all ii) £ , 

(5.7) 



(V>(ff) - </>(3, e))^( 5 )(Q e ® M)(de, d 5 ) = 0, 

i 

where \i is a measure which is equivalent to Q and $ stands for a set of functions which 
are square- integrable with respect to /x, whose linear space is dense in L 2 (R, The idea 
of the least-squares Monte Carlo method is to relax, for computational tractability, the 
principle (5.7) to 

determine w e lin ^ such that for all ib £ ^, , 
K (5-8) 

^2(ip(g k ) - (f>(g k ,e k ))ifj{g k )=0, 

k=l 

with a finite set of basis functions 

* !''..:./ 1 I\ 

and an appropriate sample 

5:=(e fc , 3fc )f =1 cM 2 , 
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chosen such that the combination -^X^fci^fefc.Sfc) °f the Dirac measures approximates 
the distribution Q £ ® fi (for instance, S being realizations of K £ N independent Q e ® fj,- 
distributed random variables). The solution to the weakened problem (5.8) is given in 
terms of 

realizations (j>(S) = (4>( e kt9k))k=i °f 4> on the sample S, 

k j ■ . ( 5 - 9 ) 

realizations M = (V^' (<?&))&= 1 j=i °^ basis functions on S, 

as follows: 

if q = (q 3 )-- =1 fulfills M T Mq = M T (f>(S), 

J 

then (5.8) is solved by ip = Qjipj- 

3 = 1 

We now formulate an algorithm for the approximate calculation of (5.5) in which the 
conditional expectation is replaced by the least-squares Monte Carlo projection. To ease 
notation, let us suppose that (£t)J = i are identically distributed (in addition to their 
independence (5.1)). 

Allowance prices via Monte Carlo method. 

1. Initialization. Given sample S = (ek: 9k)k=i ^ 2 describing the distribution of 
Q ei ® \x and a set of basis functions $ = (V'i)/=i on ^> define M as in (5.9). Set 
a T(g) — l[o,oo[(ff) f° r all g £ R and proceed in the next step with t := T — 1. 




Figure 1. The functions at for t = T— 1, . . . , T — 6, from the least-squares Monte Carlo method. 
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2. Iteration. Define a° = at and proceed in the next step with n := 0. 
(2a) Calculate </>™ +1 (S) := (a t +i (g k ~ <h « (g k )) + e fc ))f =1 . 
(2b) Determine a solution q n+1 E R J to M T Mq n+1 = M T (j) n+1 {S) . 

(2c) Define a t " +1 :=E/ =1 ^ + Vi- 

(2d) If max^Lj^ \a" +1 (gk) — af(gk)\ > £, then put n := n + 1 and continue with step 
2a). 

If max^L 1 |a" +1 (<? fc ) - a'?(g k )\ < e, then set t := t — 1. If t > 1, go to step 2, 
otherwise finish. 

Example. To illustrate allowance price calculation via the Monte Carlo method, we 
consider the following numerical example. Suppose that the penalty is set at ir = 100 and 
that the martingale increments (et)tLi are independent, identically Normally distributed. 
Note that by an appropriate choice of the emission measurement scale, the standard 
deviation can always be normalized, thus we have assumed that each e t is A/"(0.5, 1)- 
distributed. Further, consider the basis consisting of piecewise linear hut functions 

ipj-.R^- [0,1], 14 (1 - \zj -x\/h) + for xER, j=l,...,J, 

where the peaks Z\ = — (J — 1) * h/2, . . . , zj = ( J — 1) * h/2 are chosen to be equidistant 
with the distance h > 0. For numerical illustration, we set J = 16 and h := 1. Further, 
the sample S = (ek,gk)k=i f° r the Monte Carlo method is generated with K = 1000 
outcomes. For (efe)^L 1; we followed a natural choice, taking realizations of K indepen- 
dent A/"(0.5, l)-distributcd random variables. However, since the distribution of G t is not 
known in advance, an appropriate candidate for [i seems to be the uniform distribution 
concentrated on the interval which is relevant for calculations. That is, the outcomes 
(fffc)fcLi are constructed by equidistant sampling of [z%, zj], ranging from g\ = z\ = —7.5 
to gx = zj = 7.5. For the cumulative volume function c: R — » R, a H> 0.1-^/(0)+, we ob- 
served a fast and stable convergence which gave a reasonable outcome within a few 
iterations. The resulting functions (at)J=T-i are depicted in Figure 1. 

Let us outline a valuation procedure for a European call on emission allowance price. 

Valuation of European call via Monte Carlo method. 

1. Given basis functions \P = (ipj)j =1 and a sample S = {e-k,gk)k = i C IR 2 which ap- 
proximates Q £l determine (at)t = T m terms of basis coefficients using the above 
least-squares Monte Carlo algorithm. 

2. Given maturity time r E {1, . . . , T} of the European call, determine its pay-off f£ := 
(a T — K) + . Calculate least-squares projections, recursively processing for u = r, . . . , t 
as follows: 

(a) put (j)(S) = {fl{g k - c u (a u {g k )) + e k )y< =1 ; 

(b) obtain q as solution to M T Mq = Mtfi; 

(c) set =E/=i?jV ) j; 

(d) if u — 1 = t, then finish, otherwise set u : = u — 1 and return to (a). 
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3. Given recent allowance price a, calculate the state variable g as solution to a = at (g) ■ 

4. Plug in the state variable g and into function / T (i, •) to obtain the price of the 
European call as as f T {t 1 g). 

Let us conclude this section by sketching core ideas on continuous-time modeling. Our 
analysis shows that the risk- neutral allowance price evolution (At)f =0 must be described 
by a martingale whose terminal value is digital and depends on the intermediate values 
(see (3.19)). Suppose that the compliance period is given by an interval [0,T], such that 
all relevant random evolutions are described by adapted stochastic processes on 

filtered probability space (fi, J 7 , P, (J- t)te[o,T] ) , g 1Q , 

equipped with probability measure Q ~ P, 

where Q represents the spot martingale measure. Given a random variable £t and ap- 
propriate non-decreasing and continuous abatement functions Ct : R+ xf!-> R + indexed 
by t G [0,T], we follow an analogy to discrete time and consider solutions (^4t)tg[o,T] to 

A t =7r^(l {£r _ /(fct(Aa)ds > 0} |J- t ), te[0,T}. (5.11) 

Our results from the discrete-time setting suggest that if 

{the increments of the martingale (£ t = -E^^Tl-^t^telo.T] are 
independent and the abatement functions c< : R+ x — > M + (5-12) 
are deterministic and time constant (ct = c) t< £[o,T], 

then a solution to (5.11) should be expected in the functional form 

A t = a(t,G t ), te[0,T], 

with an appropriate deterministic function 

a:[0,T]xl^R, (t,g) i-> a(t,g), (5.13) 

and a state process (G() te [o.r] given by 

G t =£ t - f c s (A s )ds, f e[0,T]. 
To illustrate how such an approach allows one to guess a solution, assume that 



(O, F,P, (J"t)t G [o,T]) supports the process (W t , J r t ) te [ ,T] 
of Brownian motion with respect to Q ~ P. 

Furthermore, we respond to (5.12), supposing that 

d£t =cr t dWtwith pre-specified deterministic (o~t)t£[o,T] and 
continuous and non-decreasing abatement functions (ct = c)te\o,T}- 



(5.14) 
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To ensure the martingale property of {A t = a(t, Gt)) t6 [ 0j T] , apply the Ito formula 
dA t = da(t, G t ) = d m a(t, G t ) At + d {0>1) a(t, G t ) dG t + \d m a{t, G t ) d[G] t 
= 9(i )0 )a(t, G t ) dt - 9 (0i i)a(t, G t )c(a(t, G t )) dt + |9( , 2 )Q!(t, GtK 2 dt 
+ d {0A) a(t,G t )a t dWt 
and claim the function a as a solution on [0,T[xR to 

9(i,o)a(t,g) - d {0A) a{t,g)c(a(t,g)) + ^d {oa} a(t,g)<jf = (5.15) 
with boundary condition 

a(r iff )=7rl [0l oo[(ff) for aU 5 GM, (5.16) 

justified by the digital terminal allowance price. Having obtained a in this way, we 
construct the state process as the solution to the stochastic differential equation 

dG t =d£ t -c(a(t,G t ))dt, G =£ , (5.17) 

and then determine 

A t :=a{t,G t ), te[0,T}. (5.18) 
Finally, this process must be verified in order to solve (5.11). 

6. Conclusion 

This article explains the logical principles underlying risk-neutral modeling of emission 
certificate price evolution. We show that within a realistic situation of risk-averse market 
players, there is no connection between social optimality and market equilibrium, but 
there is a useful feedback relation characterizing risk-neutral allowance price dynamics. 
Expressing this result in terms of fixed point equations on the level of martingales, 
we address the existence of its solution and elaborate on its algorithmic tractability. 
Furthermore, we suggest an extension of these concepts to continuous time and show that 
promising results can be obtained using diffusion processes. Here, emission allowances 
and their options can be described in terms of standard partial differential equations. 
Although option pricing in this framework seems to be appealing, we believe that it is not 
superior to our Monte Carlo method since the latter can be used in high dimensions and, 
more importantly, in the presence of jumps in the martingale (£t)f = Q. This is particularly 
important to describe price shocks, which may result from possible discontinuities in the 
information flow. 
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